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Abstract 



Using a variational method it is shown that for magnetic fields B > 10 G 
there can exist a molecular ion H^ + . 
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Many years ago the formation of unusual chemical compounds in the presence of a strong 
magnetic field whose existence is impossible without strong magnetic fields was predicted 
by Kadomtsev and Kudryavtsev |I| and Ruderman [[| (for a recent advances and a review, 
see, for example, HQ and references therein). In particular, using a semi-classical analysis 
it was shown that the influence of the strong magnetic field leads to the appearance of linear 
molecules (linear chains) situated along magnetic lines. The transverse size of such systems 
should be of the order of the cyclotron radius p ~ B^ 1 ^ 2 (a.u.), while the longitudinal size 
remains of a molecular (atomic) order. These systems are called needle-like. An important 
consequence of such a quasi-one-dimensionality of Coulombic systems is the possibility of 
effectively compensating the Coulombic repulsion of nuclei. 

The goal of this Letter is to present the first quantitative study of the molecular ion + 
in a strong magnetic field which provides theoretical evidence that such a system can exist 
in magnetic field B > 10 11 G. Our study is limited to an exploration of the ground state. 
Throughout the present work it is assumed that the Born-Oppenheimer approximation holds 
which implies that the positions of protons are fixed. Exactly as for the configuration we 
consider corresponds to the case when the three protons are aligned with the magnetic field 
(linear chain, see above). Spin effects (linear Zeeman effect) are neglected. The magnetic 
field ranges from up to 10 13 G, where it is assumed that a contribution of relativistic 
corrections can still be neglected (for a discussion see, for instance, Q and references therein). 
Finally, it is also demonstrated that the molecular ion is the most bound one-electron 
molecular system in a constant magnetic field. 

The present calculation is carried out in the framework of a variational method using a 
unique simple trial function equally applicable to any value of the magnetic field strength. 
Very recently ||, this strategy was successfully applied to study the ground state of the 
molecular ion H2 and a simple 10-parameter trial function allowed one to get the best 
(lowest) values of the ground state energy for magnetic fields from up to 10 13 G (except 
B = 0, where the relative accuracy was about 10 -5 in comparison with the best calculations). 
It turned out that it was not only the best calculation in the region of accessible magnetic 
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fields — 10 13 G but the only calculation which tackled the problem using a single approach. 
Since the key point of a successful study is a wise choice of trial functions, we give a detailed 
description of how to choose trial functions adequate to the problem in hand. 

A constructive criterion for an adequate choice of trial function was formulated in || 
and further development was presented in J7|||. In its simplest form the criterion is the 
following. The trial function ^(x) should contain all symmetry properties of the problem 
in hand. If the ground state is studied, the trial function should not vanish inside the 
domain where the problem is defined. The potential Vt(x) = , for which the trial 
function is an eigenfunction, should reproduce the original potential near singularities as 
well as its asymptotic behavior. The use of this simplest possible recipe has led to a unique 
one-parameter trial function, which in particular, made it possible to carry out the first 
qualitative study of the ground state of the hydrogen molecule H 2 in the region of both 
weak and strong magnetic fields ||. Later a few-parameter trial function was proposed 
for a description of the hydrogen atom in an arbitrary magnetic field, which led, for the 
low-excited states, to an accuracy comparable with the best calculations |8].|10| • 



Now we wish to apply the above recipe to the ion H% + . Let us first introduce notation 
(see Fig.l). We consider three attractive identical centers of unit charge situated on the 
z-axis at origin and at a distance i?_ , R + from the origin, respectively. The magnetic field 
of strength B is directed along the z axis and ri 2 3 are the distances from the electron to the 
first (second, third) center, respectively. The quantity p is the distance from the electron 
to the z-axis. Through the paper the Rydberg is used as the energy unit. For the other 
quantities standard atomic units are used. The potential corresponding to the problem we 
study is given by 

2 2 2 222 B 2 p 2 

~ RZ + R~^ + R_ + R + ~ n~ 7 2 ~ 7 3 + 4 ' ^ 

where the first three terms have the meaning of the classical Coulomb energy of interaction 
of three charged centers. The recipe dictates that the trial functions should behave in 
a Coulomb-like way near the centers, correspond to two-dimensional oscillator behavior 



in the (x, y) plane at large distances and be permutationally-symmetric with respect to 
exchange of positions of the centers. It seems quite natural that the equilibrium configuration 
corresponding to minimal total energy of the system should appear at R- = R + . 

One of the simplest functions satisfying the above recipe is the Heitler-London type 
function multiplied by the lowest Landau orbital: 

\J> 1 _ e -ai(ri+r 2 +r 3 )-l3 1 Bp 2 /4^ ^) 

(cf. Eq. (2.2) in ||), where a±, f3\ are variational parameters. It has a total of four variational 
parameters if the internuclear distances R-, R + are taken as parameters. It is quite natural 
from a physical viewpoint to assume that a function of the Heitler-London type gives an 
adequate description of the system near the equilibrium position. The potential Vi(x) = 
, corresponding to this function is: 



Vt = 3al-B(3 1 -2a 1 [- + - + -\ +5^ + 2^ 

n r 2 r 3 4 



1 'p 2 + z{z + R + )\ + 



r x r 2 



— ( p 2 + z{z - R_)) + — ( p 2 + {z - R_){z + R A 
r 2 r 3 \ J rir 3 



+ a 1 p 1 Bp 2 [- + - + -] (3) 
n r 2 r 3 



It is clear that this potential reproduces the original potential (|l|) near Coulomb singularities 
as well as at large distances, \x,y\ — > oo. 

The Hund-Mulliken-type function multiplied by the lowest Landau orbital is another 
possible trial function: 

^ 2 = (e~ a2ri + e~ a2r2 + e -W^\ e -foB?/A ^ ^ 



(cf. Eq. (2.4) in [^]), where a 2 , f3 2 are variational parameters. It is obvious that this function, 
in the absence of a magnetic field, gives an essential contribution to a description of the 
region of large internuclear distances. The calculations we have performed show that this 
property remains valid for all magnetic fields up to 10 13 G. Like Eq. (0), the trial function 
(4) is characterized by four variational parameters. This function, when both internuclear 
distances are large, corresponds to a decay H 3 + — > H + p + p. 



Another trial function supposedly describes a decay mode H% + — > + p and could be 
taken to be of the form 

where c^,/^ are variational parameters. Finally, it will become obvious that the function 
(H) does give the dominant contribution to the large internuclear distances. Eq. (||) also 
depends on four variational parameters. 

To take into account both equilibrium and large distances, we use an interpolation of 
Eqs. (2), (4) and (5). There are three natural approaches to interpolate: 

(i) a total non-linear superposition: 



VTY o~ a4ri-a 5 r2~aer3 , -a 5 r\ -oaT2 — a§r-i i „— 04T\ — 0^2-0^3 i 

^A-nis-t — \e +e +e + 

e -a 6 ri-«4r 2 -a5»'3 _|_ e -asr 1 -a 6 r2-a 4 r 3 _|_ ^0^-0^2-04,^ j g-/3 4 _Bp 2 /4 ^g-j 



(cf. Eq. (2.5) in 0), where 04,5,6, At are variational parameters. The function (^|) is a 
three-center modification of the Guillemin-Zener type function used for the description 
of the molecular ion if 2" in a magnetic field ||. If = = = a±, the function 
(H) reduces to Eq. (2). When = ct2, a$ = = 0, it coincides with Eq. (4). Finally, 
if a 4 = a 5 = «3,a 6 = 0, the function @ reduces to Eq. (§). In total there are 6 
variational parameters characterizing the trial function (H); 

(ii) a partial non-linear superposition: 

this appears if in Eq. (||) the two parameters are equal, for instance, a 4 = a 5 : 

iT; _ [ „-04(n+r 2 )— a 6 r3 1 „-04(n+r 3 )— a 6 rz , ^-cuC^+^-aen \ -,&4-Bp 2 /4 

^4-nls-p — l e ~r e ~r e l e 5 

(7) 

This function can be considered as a non-linear interpolation between Eqs. (4) and 
(5). 
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(iii) a linear superposition of Eqs. (2), (4), (5) 

V 4 _ ls = A^i + A 2 ^ 2 + A 3 ^ 3 , (8) 

where the relative weights of Eqs. (2), (4), (5) in Eq. (7) are taken as extra variational 
parameters. This is a 10-parameter trial function. 

Of course, as a natural continuation of the above interpolation procedure one can take a 
linear superposition of all five functions (2), (4), (5), (6), (7): 

—nls—t 4—nls—t 

= M-nls-t^i-nls-t + A 4 _„^- p ^ 4 -n/ S - P + ^1^1 + A 2 ^ 2 + A 3 ^ 3 , (9) 

(cf. Eq. (2.7) in ||), where again, as in the case of the function (8) the relative weights 
of different, 'primary' trial functions are considered as variational parameters. In total, the 
trial function (9) is characterized by 17 variational parameters. However, only part of our 
calculation is carried out using this function. Usually, some particular cases of Eq. (9) 
are explored. The general case will be presented elsewhere. The minimization procedure 
is carried out using the standard minimization package MINUIT from CERN-LIB on a 
Pentium-Pro PC. All integrals were calculated using the CERN-LIB routine DGAUSS with 
relative accuracy < 10~ 7 . 

In Table I the results of our variational calculations are presented. It is quite remark- 
able that for magnetic field strengths > 10 11 G there exists a minimum of total energy in 
the plane. Furthermore, for such magnetic fields the value of the energy at the 

minimum correspondent to the total energy of H 3 + is always lower than the total energy of 
the hydrogen atom but higher than that of H 2 . Hence the decay mode H% + — * H + p + p 
is forbidden. However, -£^3~ + is unstable with respect to the decay H 3 + — > H 2 + p. It 
seems natural to assume that even if one-electron systems like Hf ++ , H§ +++ etc. would be 
bounded, their total energies will be larger than the total energy for H 3 + . This assumption 
and comparison of the total energies of H, H 2 , H 3 + (see Table I) allows one to conclude that 
H 2 is the most stable one-electron system in a constant magnetic field oriented along the 
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magnetic field. The equilibrium distances for H£ + decrease with the growth of the magnetic 
field: the ion H^ + , like H2, becomes more and more compact. It is worth noting that for 
both H2 and H% + the average value (z) is much smaller than a 'natural' size of a system 
determined by the positions of the centers: R eq for and (R +eq + R- eq ) for (see Fig. 
1). In other words this means that the localization length of electron is much smaller than 
the 'natural' size of the system. 

Fig. 2 shows the electronic density distribution as a function of magnetic field. For a 
magnetic field B ~ 10 11 G the distribution has three clear maxima corresponding to the 
positions of the centers, but the electron is situated preferably near the central proton. The 
situation changes drastically with an increase of magnetic field: the electron is localized 
near z — 0, having almost no memory of the two centers on either side. It is important 
to investigate paths of possible tunneling. There are two explicitly pronounced (symmetric 
with respect to R + <-> R_) valleys in the electronic potential energy surfaces, E tota i(R + , R_) 
(see, Fig. 3a). 

They vary from the position of the H$ + minimum to infinity which corresponds to the 
H£ +p dissociation: (R +eq ,R_ eq ) -> (oo,R e £+) and (R +eq ,R- eq ) -> (R e £+,oo), where R e £+ 
is the equilibrium distance for the ion. In Fig. 3b one can see the profile of the valley 
as a function of magnetic field. Calculating the Gauss curvatures in the H% + minimum 
one can estimate the position of the ground state energy level and answer the question of 
whether the well is deep enough to hold an energy level. It is always delicate to answer this 
question starting from what 'depth' of the well the level exists. Usually, it requires special 
analysis. We made an estimate and obtained the result that for a magnetic field of 10 11 G, 
the situation is not certain, the well is probably still too shallow to hold the ground state 
energy level. However, the well undoubtedly becomes sufficiently deep for 10 12-13 G. From 
the form of the profile (see Fig. 3b) it is quite obvious that for 10 12 ~ 13 G the barrier is rather 
high and the probability of tunneling should be small. 
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TABLES 



B = 10 11 G B = 10 12 G 

Leg E R eq (z) E R eq (z) ^ ±v eq 

(Ry) (a.u.) (Ry) (a.u.) (a.u.) (Ry) (a.u.) (a.u.) (Ry) (a.u.) (a.u.) 



B = 
E R P 



B = 10 13 G 
E R P 



H -1.000 36.929 413.57 4231.6 

H+ -1.205 1.997 35.036 0.593 0.312 408.300 0.283 0.174 4218.662 0.147 0.107 

H+ + 36.429 0.803 0.432 410.296 0.346 0.219 4220.090 0.165 0.121 



TABLE I. Data for the ground state of and a comparison with data for other one-electron 
systems, H, ■ Total energy E is in Rydbergs, the equilibrium distance R eq = R+ eq = R- eq 
(see text) and the average value of the longitudinal size of the system (z) in a.u. Total energy for 
hydrogen atom from [11]; data for H2 from [5]. 
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FIG. 1. + in a magnetic field B. Explanation of the notation used 
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FIG. 2. Electronic distribution for various magnetic fields: 10 11 G (a), 10 12 G (b) and 10 13 G 
(c). It peaks more and more sharply at origin with growth of magnetic field. 
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FIG. 3. Valleys in the electronic potential energy surfaces, E tota i(R + , R-), (a) and their profiles 
(b) for various magnetic fields. The position of the minimum is indicated by a bullet, and AE 
denotes the depth of the well: the distance between top of the barrier and the value of minimum. 
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